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Чисельний анал!з конвективно! модел! кристал1заци 


Исследуется пространственная задача Стефана с учетом примеси и конвективного движения в жидкой фазе, 
описываемого уравнениями Навье-Стокса. Предложен метод изучения этой задачи, состоящий в разложении 
решения в ряд по степеням малого параметра. В нестационарном случае решение соответствующих краевых 
задач для определения членов разложения строится как неподвижные точки операторов. Исследовано 
влияние конвекции на фронт кристаллизации. Разработан также метод решения задач сопряжения, 
возникающих при исследовании задач Стефана в пространстве. 

Ключевые слова: дифференциальное уравнение, свободная граница, численные методы, 
функционал, оптимизация. 


Тре @гее-4итепз1опа| сопуесноп 5{еРап ргоМет ш Па рВазе 1$ шуезисжеа. ТЬ$ рго ет 1$ дезси6еа 
Бу Мамег-ЗюКез едааНопз. ТВе тео Гог гезеагсЬ оЁ 1$ рго ет, мс соп$15$ оЁ фе зоаноп 
ехрапз1оп Гог зетез оЁ а зтаП рагатеег, 15 оНегед. ш поп- %“айопагу сазе, Ше 4ес1$1оп оЁ соггезропд те 
Боипдагу-уа[ае (азКз Гог дейшиоп оЁ Ше ехрапзюп тшетБегз 15 ЮюгтеЯ аз Нхе ро оЁ орегаогз. 
Сопуесйоп шЙчепсе оп Ше Шоп оЁ сгубаШтайоп 15 ехр]огед. Тве тефо4 ог зоаНоп оЁ Фе сопиганоп 
(азКз, ушсВ паке р1асе \Веп зва4уте ЗеРап рго ет 1т зрасе, 1$ а1зо деуеюореч. 

Кеу мога$: аШегепна] едиайоп, ее Боипаагу, питегса] а]2огИтп$, Капс@опа|, орипихайоп 


Дослджуеться просторова задача Стефана з урахуванням домппок 1 конвективного руху в р1диннйй фаз, 
як! описуються рвнянням Нав’е-Стокса. Запропонован метод вивчення ще! задачу, який полягае в розклад! 
розв’язку у ряд, ыдповдно до ступенв малого параметру. У нестацюнарному випадку розв’язання 
в1дпов!дних крайових задач для визначення членв розкладання будуються як нерухом! точки операторв. 
Дослджено вплив конвекцй на фронт! кристалзаци. Розроблено метод розв’язання задач спряження, що 
виникають при дослдженн! задач Стефана в простор. 

Ключов! слова: диференщальне р1вняння, в1льна межа, чисельн! методи, функщонал, оптим1заця 


1. Статья посвящена численной реализации нелинейной конвективной задачи 
теплообмена, возникающей при кристаллизации вещества, в следующей постановке. 


Пусть О, — заданная область в К°, граница которой состоит из двух замкнутых 


связанных гладких поверхностей Гу и Г., не имеющих самопересечений. Пусть, далее 
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Г, - гладкая замкнутая поверхность, лежащая внутри ©), ‚ такая, что Г, лежит внутри 
ограниченной области, границей которой является Г, . Поверхность Г. разбивает 
О, на две подобласти (№ и О’, которые в начальный момент { = 0 заняты жидкой и 


твердой фазами соответственно. Будем обозначать через О" область занятую жидкой 
(твердой) фазой в момент времени 1. Заметим, что в процессе кристаллизации проходит 
изменение границы Г’, (это связано с тем, что жидкая и твердая фазы имеют разные 


плотности), а граница Г. остается неизменной. Задача состоит в определении 
областей © и ©. (т.е. границ Г’ и Г,), занимаемых твердой и жидкой фазами 


соответственно в момент времени 1е[0,Т], вектора скорости У(х, г) = 


=(У(х, 1), У, 1, У, (х,1)), давления р(х,г), концентрации примеси с(х,!), распределений 


температур жидкой и’ (х,ё) и твердой и (х,#) фаз по следующим условиям: 


Як +(УУ) и (ва Уи (к, =0, (х,1) Е 0%; ны: —а\и (х0=0, (х,) ЕР; 
9 ›1) +(УУ) У(х, УР Г) =У\У?У (С + " (и’,с), УУ(х, Юй =0, (ЕР 
У(х,0) = ССЭ; ТУ, ри =, (хр) Е Г: У, =а- ем, У =0, (х,0 ЕГ, 
ги =вои боеггог ти = А*(х); (1) 
и*=и` =Т'-5С,К. А т ро еГь 
дп [61 
Ос(х, г) 


о ^^ + (УУ)с(ь, й- —УУ°с(х, 0), (х,)ЕДр’; с(х,0) = 2. (х), 


с(х, г) = 8(х,1), (х,) Е Г’; о - Всй’,, (х,)ЕГ,. 
и 


здесь 0»; ={(х,0): хе О», Ее (0,Т)}, х=(х,х,,х,), О* — области соответственно жидкой 


и твердой фаз, д0*=Г,оГ*, 00-=Гуог,, п - нормаль к Г, , направленная в 

сторону ©*; Ту, в, Хх, р’, ра, В, 1, к,к,, - Положительные постоян- 

ные, у = (-°_ о - © ›, Т(У, р) — тензор напряжений с элементами т, = оо у 
Ох, ’дх, ’дх, Хх, *бх 


У’, — скорость движения фронта кристаллизации в направлении нормали п, У иу - 


нормальная и тангенциальная составляющие у. 


Если Ф(хр=и’(хЮ-+Ес(хр)-Т’=0 - уравнение поверхности Г,, тогда 
и, =-Ф,/|УФ |, р ВА 
[У (иг +а=с) | 


Укажем, что условие Стефана можно представить также в виде: 
к? |Уи  [^-к? [Уи* [+= (к? +К К, (Уи ,Ус)-=(к? + КК, )(Уи*, Ус) + хр (к_и; +к. и’) + 
+хр'Е (к, +к_)с, =0, (х,Г)е Г,. 


2+9, 
Предполагается, что АСЕ Н““(О”), СЕН” (О), Во) ен“? Гат), 
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2-+а, 
(и, се СКВ”), в(хеЕН “? (Г*х1О,Т]), 8) ЕН““ (5). При этом 2(х,1) и 
1+4 


14а,—— 
8, (х,Р) должны быть функциями класса Н `? (Е®*х[О,Т]). Предполагается также 


выполненными условия согласования до первого порядка включительно, которые 
следуют из предположения существования гладкого решения и формулируются 
аналогично [1 с. 363, с. 268]. 

Отметим, что при малых значениях [, задача (1) разрешима в классе гладких 


+ ой Е — ыы =. ры ЕЛЕ в ЕЕ 
функций и ЕН *(О;), УЕН ? (Р:), СЕН *(2;), УреН *(О;), а 


2-9 


2+4 , 
границы Г’ и Г, описываются функциями, принадлежащими классам Н ®. [2]: 


Решение задачи (1) моделирует процесс кристаллизации вещества с учетом 
переноса примеси в жидкой фазе. При этом последнее условие в (1) следует из 


закона Нернста, а }(и’,с) описывает влияние неравномерного распределения тем- 
пературы и концентрации примеси на движение жидкости. 
2. Известно, что свободные границы Г, и Г’ можно представить в виде 
5 у * * = * * * * 
Г; ={х = х(®) +п(®)р(©,1)}, Г’ ={х=х(® )+1@® по )}, где в =(@,0,), ® = (®, ,®,), 


х(®о) Е Гу, х(®’)е Гу, р(®,В и 1 (©’,Г) некоторые функции соответственно классов 


2+ 2-а, 


2-0, ыы 2+ —— % 
Н 2? (Г. хо] и Н 2 (Го х[0,Т]), р(©,0) =Оип(® ,0)=0 [2]. 
Предложен метод решения задачи (1), состоящий в разложении решения в ряд 
по степеням малых чисел 5: 


и (х,1:5) = (0+ ='и(х,) ‚ рот; &) = рсд+У =" р, (х,Г), 


= = 


У(х, 58) =, +У =, (х,0,1=1,2,3; р.б) =У'=“р, 0, — (2) 
К=1 


К=1 


с(х, в) = с, (%) + У: г (х,в. 


Для нулевого приближения и У. = (И 09,0, 09), Г и с (*х) из 


условий (1) и разложения (2) вытекает следующая задача: 
(ИУ, (>) +Ур, (0) =УУ?У, (2х) + Т (и ,с,), хеОй, 


УИ) =0, хеО,", ТИ, ри =-а(®я, хеГо', У, =-а-РОМ,, У =0, 
р 


хеГ,; (УУ)и* —а Уи! =0, хЕО*, иЁ (Хх) = В*(х), хе ГГ, щ(х=и (С =Т’, (3) 


ан 0“ А 
хеГо К.Е, - =0, хеГ,, У?из =0, хеО, (ИМ)с, -УУс, =0, 
+ + су 
хЕО,, с (х) = в, (х), хЕГ`; а. хеЕГ.. 
И, 
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Здесь ради простоты предполагается, что функции В* и 4 зависят только от 
переменной х. 

Лемма 1. Пусть функции и,’ (х) = А*(х), У = Со), с (х)= 8, (х) являются 
решением задачи (3) соответственно в области © и ©». Тогда эти функции можно 
взять в качестве нулевого приближения задачи (1). 

3.Далее; пусть: О =, МО Гетеро], 


7 + 
Рассмотрим первое приближение (У,и,,р,р.,с,) задачи (1) для малых чисел 


5. Имеем: 
ду, : й ” 2 217 и, + о ЗИ + 
5, + (ИУ)У, + (ИМИ, + Ур, =УУТУ, + Х, (ши, + Кс (му, сис, (х,ПОЬ; 
УУ =0, бьрЕО- (4) 
Т(У, +У,р)п=0, хЕГ,*,У (х,0) = 0, 
ау О жеОа 
р |Ущ | 
и! (уу) + уу в 272 + .. ди: 2 2. - _ _ 
- 1У)мо + (ИУ)и, - а. Ум, =0,(хЕО.; 5: +а\Ути, =0,(х0еЕО., 
и* (х,0) = О; и (х,р)=0, (хр ЕГу Г, (5) 
ь - ди. ди* ао 
и — Е +) = Хр ет О ЕГи; 
Ос, ту. Ту” 2 + 
О ОНО 
с. (х,0) = 0, с, (х,6) =0, (хе Гуг, 
д - (6) 
С и, 
= 4х, ое ет: 
дп | Ум, | 
АЕ (.В=0, (хоеГе. 
п 


Зададим теперь У= У (хп. Затем решим задачу (5), (6) и найдем и’,с,р. После 
чего заменим и’,с,р — решением задачи (5), (6) и решим задачу (4), являющуюся на- 
чально-краевой задачей для системы Навье-Стокса. Затем, используя новое значение 
У(х,г), снова решаем задачу (5) и (6) и тд. Таким образом, получим процесс последо- 


вательных приближений. Доказательство сходимости этого процесса аналогично приведен- 
2+9 ый 
ному в работе [3]. При этом при заданном р, (®,еН * (Гг) найдем функции 
2+9 2+9 


НЕЕ 2+9,—— Е 
и(х,5р)ЕН * (0), с (хЕр)еН ' (0.), как единственное решение задачи 
(5), (6) [1], причем р, (©,Г) находим как неподвижную точку сжимающегося оператора 
ди: й ди’ 
дп ‘дп 


М, : М р, = — [| + /(х. 6), ХФ) еГу,. Имеют место следующие утвер- 
0 


ждения. 
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д 
Лемма 2. Пусть выполнено условие | УА* (х) |= чи на Г, . Тогда оператор М, , 
и 


2+9 2+9 
2+а,—— 2+9, —— 
действующий из Н °(Г;)вН —? (Гу), имеет там неподвижную точку. 


Лемма 3. В качестве первого приближения задачи (1) можно взять решение за- 
дачи (4) — (6): и*(х,в),с, (ХВ, (х,8, р, (ХВ, р, (%0. 


08, (х 
Теорема. Пусть Нм, О на Гу’. Тогда при малых числах & и достаточно 
т 


малых значениях { справедливы формулы: 


№ и: (х(6), 1) 


Г ее. (хрЕГ 


с (х(6),1) + 8, (х(®))- 8(х(6),1) 
08,(х(6)) 
Оп 


Г*:х=х(®’)-=п +0(=), (х,#) Е Г* 


07 › 


2+9 
+а,—— 


2 
где и: (х,1), с, (х,0), р, (®,0), 1, (©,1) — функции класса Н ? , являющиеся 
решением задачи (4) - (6). 


4. Рассмотрим случай, когда © = в < в} и без учета кон- 


центрации примеси в задаче (1). Далее, так как решение задачи (1) ищем в виде (2), 
тогда нулевое приближение находим как решение следующей задачи: 


Я = В* (х),и* (х) |, =0, ( 
Уи (х)]- Ум (х)]=0, хе Г, 7 


Уи" (х) =0, хе Ох, А+ (х) 


С(х)=0,хеО:; 


Заметим, что замена й =Ки, при хеО ий’ =и', если хеО`, сводит задачу 

(7) к случаю К =1. Поэтому в дальнейшем будем считать, что это условие выпол- 
+ ` 

нено. Нулевое приближение и’(х), Г, найдем из условия минимума функционала 


(ие. Г) = [Ум Фак 4хь, здесь @ = 0О. и и=и при хЕО ии=и,, если 
[Ф) 


же”. 


Далее, рассматривая функционал Г в сферических координатах, получим 
ы тт К й 1 ь 1 ь г 
Ци) = | [а > а } т @афа®ар. 


Минимум функционала ищем в следующем виде 


ив (в--В*)+(* р?) (р?- п) Ср". 
В к к=0 . 


Неизвестные коэффициенты С определяются методом Ритца. В частности в 


случае нулевого приближения 
2 2 


о ВР (в--В")+(В*-р?)(р? 72), 
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из уравнения ОГ (м)/ 0С, =0 определим коэффициент С... Свободную поверхность 
Го: Р=Ро (ф,9) найдем из условия и, (0, Ро (х,9)) =0. Тогда для свободной по- 


верхности /’, получим уравнение 


и; (<, 9,1) 


Г, :р(Ф,©9,1) = р. (ф,0)-Ке— 
( ) о ( ) УА* (Ф.0) 


+0(Ве). 


На рисунке представлена поверхность Г ‚ При следующих значениях параметров: 


#=200, К=6, г=0,8, -^<0<7, -Т 557, 
2 3 2 р) 


2 2 — 2 2 
В* =3| с0$ 9+с0$ Ф |, В =-0,35] сз 9 -с0$ Ф |--0,1 
Предложенный алгоритм построения поверхности Г р позволяет исследовать 
эту поверхность в зависимости от основных параметров задачи (1), в том числе и 
числа Ке. На рис. 1 свободная поверхность Г, расположена между сферами ра- 


диусов А И г. 


Рисунок 1 
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А.5. МтепКо, 5.А. СупКко 
№Митептса[! Апа[у515 Моае о7 СтуяаШсаноп тий Сопуеспоп 


ЗП, ошу опе рарег Ваз Бееп Кпо\уп аз фай опе, \умсВ сопатз$ Фе ргооЁ оЁ Ше ех$епсе оЁ Ше 
с1а5$1са[ зо[аноп ю Ше м№о-рВазе тшШа-дптепз1опа| $(еГап ргоет: 


ы — 

+ буди р) - а? У?и* (хр) =0, (хе, 
т аЗу?и (д =0 (хрЕРх, 
РР + буи, + Урон = ПУ, + и) 

ё 


и" (х, Е) 


нар =0, >. к. о К. ий со5( и, х,) + К с0$( и,1) =0› 
и Ох, дх, 


+ 
+ + — 
рт = {хх ЕО; ‚1 > (0,2) х = (х1,х2,хз), О, . 
Тре олуеп агафе 1$ Фе зесоп4 опе, лев сопаш$ Ше гезай оф по 1е5$ 1етийсапсе. ТВе аифог$ 
сопза4егз Ше поп-4есепегае г\о-рВазе 5{еРап ргоет т а суйпаег, аззитлте фа Феге ех15$ а зооп оЁ 


+ 
Ше соггезропате зайопагу Зап ргоМет И() (х) ап Ша Фе Боипдагу ап Фе пийа| даа оЁ Фе 


ргоет аге ошу $П=ВИу ЧШетгеп( бот Фе соггезропа те свагасет$с$ оЁ Фе ханопагу зоайоп: 


и= (х,1;Ве) = иё (х) + у (ве) бы) 
К =1 


о . 
У; (х, Ве) =Уц(х) + У (Ве) У (х,Р), в =1,2,3; 
К =1 


р (®,1;Ве) = У (Ве “рев ,0.. 
К =1 


Тре ехт$епсе оЁ Фе зтоо@ зоайоп оЁ Ше ргоШет 15 ргоуе4 ипаег по юо гезблсйуе аззитриоп$ 
сопсегише Фе шри аа оЁ Фе ргоет. Могеоуег, Фе аифогз$ ш@сме а тре соп@ оп, Усв 
сиагащеез Це ех1%епсе оЁ Фе соггезроп4 те зайопагу зойоп апа из забИиу, Ше ]1айег Бет ип4егзюо4 
Ш 15 изиа[ зепзе. ТЫ$ рарег тау Бе зиоп?]у гесоттепае4 №ю аП пиегежме4 шт рагабойс Нее Боипдагу 
ргоетл$ ш хепега| ап4 т Фет с1азз1са| теайтепе ш рагасЧаг. 


Статья поступила в редакцию 19.12.2011. 
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